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Abstract — We define a class of multi-hop erasure networks that 
approximates a wireless multi-hop network. The network carries 
unicast flows for multiple users, and each information packet 
within a flow is required to be decoded at the flow destination 
within a specified delay deadline. The allocation of coding rates 
amongst flows/users is constrained by network capacity. We 
propose a proportional fair transmission scheme that maximises 
the sum utility of flow throughputs. This is achieved by jointly 
optimising the packet coding rates and the allocation of bits of 
coded packets across transmission slots. 

Index Terms — Code rate selection, cross layer optimisation, 
network utility maximisation, packet erasure channels, schedul- 
ing 

I. Introduction 

In a communication network, the network capacity is shared 
by a set of flows. There is a contention for resources among 
the flows, which leads to many interesting problems. One such 
problem, is how to allocate the resources optimally across 
the (competing) flows, when the physical layer is erroneous. 
Specifically, schedule/transmit time for a flow is a resource 
that has to be optimally allocated among the competing 
flows. In this work, we pose a network utility maximisation 
problem subject to scheduling constraints that solve a resource 
allocation problem. In another work, we studied the problem 
of optimal resource allocation in networks [1]. 

We define a class of multi-hop erasure networks, and 
consider packet communication over this class. The network 
consists of a set of C > 1 cells C = {1, 2, • • • , C} which 
define the "interference domains" in the network. We allow 
intra-cell interference {i.e transmissions by nodes within the 
same cell interfere) but assume that there is no inter-cell 
interference. This captures, for example, common network 
architectures where nodes within a given cell use the same 
radio channel while neighbouring cells using orthogonal radio 
channels. Within each cell, any two nodes are within the 
decoding range of each other, and hence, can communicate 
with each other. The cells are interconnected using multi- 
radio bridging nodes to create a multi-hop wireless network. 
A multi-radio bridging node i connecting the set of cells 
B{i) = {ci, .., Cn} C C can be thought of as a set of n single 
radio nodes, one in each cell, interconnected by a high-speed, 
loss-free wired backplane (see Figure [T]). 

This work is supported by Science Foundation Ireland under Grant No. 
07/IN. 1/1901. 
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Fig. 1. An illustration of a wireless mesh network with 4 cells. Cells a, 
6, c, and d use orthogonal channels CHi, CH2, CH3, and CH4 respectively. 
Nodes 3, 5, and 6 are bridge nodes. The bridge node 3 (resp. 5 and 6) is 
provided a time slice of each of the channels CHi & CH2 (resp. CH2 & 
CH4 for node 5 and CH2& CH3& CH4 for node 6). Three flows /i,/2, 
and /s are considered. In this example, Cf^ = {a, 6}, Cf^ = {d, b, a}, and 
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{c,d}. 



Data is transmitted across this multi-hop network as a set J^ 
= {1,2,--- ,F}, F > 1 of unicast flows. The route of each 
flow / G J^ is given by Cf = {ci(/), C2(/), • • ' ,Q/(/)}, 
where the source node s{f) G ci(/) and the destination node 
d{f) e C£j^{f). We assume loop-free flows (i.e., no two cells 
in Cf are same). Figures [T] and [2] illustrate this network setup. 
A scheduler assigns a time slice of duration Tf^c > time 
units to each flow / that flows through cell c, subject to the 
constraint that ^ n^^^ Tf^c < ^c where Tc is the period of the 
schedule in cell c. We consider a periodic scheduling strategy 
(see Figure O in which, in each cell c, service is given to the 
flows in a round robin fashion, and that each flow / in cell c 
gets a time slice of Tf^c units in every schedule. 

The scheduled transmit times for flow / in source cell ci (/) 
define time slots for flow /. We assume that a new information 
packet arrives in each time slot, which allows us to simplify 
the analysis by ignoring queueing. Information packets of each 
flow / at the source node S{f) consist of a block of kf 



symbols. Each packet of flow / is encoded into codewords 
of length Uf = kf/vf symbols, with coding rate < r/ < 1. 
The code employed for encoding is discussed in Section [III 
We require sufficient transmit times at each cell along route C/ 
to allow Uf coded symbols to be transmitted in every schedule 
period. Hence there is no queueing at the cells along the route 
of a flow. It is not apparent at this point whether it is optimal 
for flow / to transmit a single code-word of n/ symbols or 
transmit a block of n/ symbols where each block carries some 
portions of each of a set of coded packets. 

Channel Model: The channel in cell c for flow / is 
considered to be a packet erasure channel with the probability 
of packet erasure being /3/,c ^ [O^l]- Thus, the end-to- 
end channel for flow / is a packet erasure channel with the 
probability of packet erasure being 

/^/=i- n [^-pfA 

ceCf 

Let the Bernoulli random variable Ef [i] indicate the end-to- 
end erasure seen by the i\h block of flow / (independent of the 
erasure seen by other blocks) of flow /. Note that Ef[i] = 1 
means that the ith block is erased, and Ef [i] = means that 
the ith block is received successfully. Note that P{Ef[i] = 
l} = /3/ = l-P{^/[i] = 0}. 

Each packet has a deadline of Df slots, by which time 
it must be decoded. Such a delay constraint is natural in 
applications such as video streaming. A packet is in error if the 
destination fails to decode the packet by the deadline. Letting 
ef{rf) denote the error probability that a packet fails to be de- 
coded before its deadline, the expected number of information 
symbols successfully received is Sf{rf) = kf{l — ef{rf)). 
Other things being equal, we expect that decreasing r/ (i.e., 
increasing the number of coded symbols Uf = kf/rf sent) 
decreases error probability e/ and so increases Sf. However, 
since the network capacity is limited, and is shared by multiple 
flows, increasing the coded packet size n/^ of flow /i gen- 
erally requires decreasing the packet size Uf^ for some other 
flow /2. That is, increasing Sf^ comes at the cost of decreasing 
Sf^ . We are interested in understanding this trade-off, and in 
analysing the optimal fair allocation of coding rates amongst 
users/flows. 

Our main contribution is the analysis of fairness in the 
allocation of coding rates between users/flows competing for 
limited network capacity. In particular, we adopt a utility- 
fair framework, and propose a scheme for obtaining the 
proportional fair allocation of coding rates, i.e. the allocation 
of coding rates that maximises '^f^jr^ogSf{rf) subject to 
network capacity constraints. This problem, which we show 
in Section III, requires solving a non-convex optimisation 
problem. Specifically, at the physical layer, the (channel) 
coding rate of a flow can be lowered (to alleviate its channel 
errors) only at the expense of increasing the coding rates of 
other flows. Also, at the network layer, the length of schedules 
of each flow should be chosen in such a way that it maximises 
the network utility. Interestingly, we show in our problem 
formulation that the coding rate and the scheduling are tightly 
coupled. Also, we show that for a log (network) utility function 
(which typically gives proportional fair allocation of resources) 




Fig. 2. An illustration of transmission scheme in cell a of the network 
in Figure [l] Every transmission schedule of Ta time units is time-shared by 
nodes 1 and 3. Note that (f)/\(f)NfRf symbols of the encoded packet p are 
transmitted in transmission schedule p-\-A, where A G {0, 1, 2, • • • , nf — 1}. 
The scheduling or capacity constraint of cell a may not be tight (indicated by 
empty time slice in the figure), as the rates of flows /i and /2 are governed 
by the whole network. 



the optimum rate allocation (in general) gives unequal air- 
times which is quite different from the previously known 
result of proportional fair allocation being the same as that of 
equal aii^time allocation ([2|). This problem, which we show 
in Section Hill requires solving a non-convex optimisation 
problem. Our work differs from the previous work on network 
utility maximisation (see fSl and the references therein) in 
the following manner. To the best of our knowledge, this 
is the first work that computes the optimal coding rate for 
a given scheduling (or capacity) constraints in the utility- 
optimal framework. 

The rest of the paper is organised as follows. In Section [III 
we obtain a measure for the end-to-end packet erasure, and 
describe the throughput of the network. We then formulate a 
network utility maximisation problem subject to constraints on 
the transmission schedule lengths. In Section III, we obtain the 
optimum transmission strategy and the optimum packet-level 
coding rates for each flow in the network. In Section |Vl we 
provide some simple examples to illustrate our results. Due to 
lack of space, the proofs of various Lemmas are omitted. 

II. Problem Formulation 

The encoding has two stages. The first stage is the encoding 
of each information packet using a standard generator matrix 
such as a Reed-Solomon code or a fountain code IJH . Let Pf [t] 
denote the information packet that arrives at the source of 
flow / in slot t. A packet Pf [t] of flow / has kf symbols, the 
encoded packet Cf [t] of which is of size n/ = kf/rf with < 
r/ < 1, and we assume that the code is such that the packet 
Pf [t] can be reconstructed from any of its kf encoded symbols 
(this is possible, for example, by Reed-Solomon codes). 

The second stage allocates the content of the encoded packet 
Ct of the first stage across the transmitted packets. Each 
encoded packet is segmented into Df portions (where we 
recall that Df is the decoding deadline requirement for each 
packet of flow /), the size of the Ath portion being (/)/(A)n/, 
where A G {0, 1, • • • , L>/ - 1} and ^ (/)/(A) ^ 1. At 
transmission slot t, a transmitted packet is assembled from the 
0/(0) portion of Cf[t], the ^/(l) portion of (7/ [t — 1], and so 
on until the (t)f{Df — l)th portion of packet Cf[t — Df + 1]. 
This procedure is illustrated in Figure [5] for n/ = 3. Note 
that the transmitted packet is of size n/ symbols. To decode a 
packet Pf [t] of flow /, we use the transmitted packets that are 
received during the transmission slots t,t + l,--- ,t + Df — 1. 
Note that the conventional strategy of transmitting an encoded 
packet every transmission slot corresponds to the special case: 
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III. Network Utility Maximisation 

For flow /, the total expected throughput as a result of 
transmitting T > 1 packets is given by 
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Fig. 3. Two stage encoding (example of Df =3): information packet Pf [1] P {<^/ [t] = X^ , t = 1 , 2 , • • • , T} is not a product-form 

of size fc/ is encoded to Cf [1] of size Uf = kf /rf , the contents of which are distribution as the packet erasures a/ [t] s are correlated, 

allocated across subpackets ^oCf[l],^iCf[l],^2Cf[l] across 3 timeslots. However, the above expectation can be written as 



^f{0) = 1 and ^f{l) = ^f{2) = ••• = ^f{Df - 1) = 0. 
We call the transmission scheme outlined above with general 
^.(A)s a generalised block transmission scheme. 

A. Network Constraints on Coding Rate 

Let Wf^c be the PHY rate of transmission of flow / in cell c. 
For each transmitted packet of flow /, each cell c e Cj along 
its route must allocate at least ^^ units of time to transmit 
the packet (or encoded block). Let J^^:={/gJ^:cgC/} 
be the set of flows that are routed through cell c. We recall 
that the transmissions in any cell c are scheduled in a TDMA 
fashion, and hence, the total time required for transmitting 
packets for all flows in cell c is given by ^f^jr^ ^^. Since, 
for cell c, the transmission schedule interval is Tc units of time, 
the coding rates r/ must satisfy the schedulability constraint 
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B. Error Probability - Upper bound 

Lemma 1. The end-to-end probability of a packet erasure 
for flow f is bounded by 



e/ 



Df-l 



kf 



^{ E M^)^Ef[A]>nf-kf 
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< exp 



=: e/. 



Df-l 



Of{l-rf)- Y. ln(l-/3j+/3je 
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where Of > is the Chernoff-bound parameter 

Let the random variable af [t] indicate whether packet Pf [t] 
is successfully decoded or not, i.e., 

af[t] 



1 0, c 



if packet Pf [t] is decoded successfully 
otherwise. 



We note here that the decoding errors for the successive 
packets are correlated, as each encoded packet overlaps with 
the transmission of previous Df — 1 packets and the successive 
Df — 1 packets. Hence, the sequence of random variables 
a/[l], a/[2], a/[3], • • • are correlated. But, the probability of 
any af [t] = is upper bounded by Lemma [T] 



SfiT) = Yl E kfXtP{aAt]=xt} 

= T-kf{l-ef) 
Thus, the (average expected) flow throughput is defined as 



lim 



Sf{T) 



= kf{l-ef). 

We are interested in maximising the utility of the network 
which is defined as the sum utility of flow throughputs. 
We consider the log of throughput as the candidate for the 
utility function being motivated by the desirable properties 
like proportional fairness that it possesses. 

We define the following notations: the Chernoff-bound 
parameters 6 := [Of]feT, coding rates r := [r/J/^j^, and 
the allocation of coded bits across transmission slots $ := 
[(t>f]feT where 0/ := [(^/(O), 0/(1), • • • ,0/(1^/ - 1)]. Thus, 
we define the network utility as 

U{^,0,r) := Y\n{kf{l-ef{cl>f,Of,rf))) 
feT 

=: ^ln(fe/) + /7(^,^,r) (1) 

The problem is to obtain the optimum coded bit allocation $*, 
the optimum Chernoff-bound parameter ^*, and the optimum 
coding rate r* that maximises the network utility. Since, kf, 
I the size of information packets of each flow / is given, 
'maximising the network utility is equivalent to maximising 
/7(^, 6^ r) := ^f^j: In (1 — e/). Thus, we define the follow- 
ing problem 
PI: 



max Ui^.O.r) 
kf 



subject to 2_\ 



f:ceCf 
Df-l 



rfWf,c 



J2^f(^y 



< Tc, \/ceC 



V/G^ 



(2) 



(3) 



4>f{A) > 0, 


V/ € J', < A < D/ - 1 


Of >0, 


V/€^ 


rf<Xf 


v/e^ 


rs>\f 


v/€ J- 



We note that the Eqn. Q enforces the network capacity (or 
the network schedulabiHty) constraint. The objective function 
[/(^,^,r) is separable in (0/, 6>/, r/) for each flow /. 
Importantly, the component of utility function for each flow 
/ given by ln(l — 6/(0/, 6^/, r/)) is not jointly concave in 
((/)/,6>/,r/). However, In (1 — 6/(0/, 6>/,r/)) is concave in 
each of 0/(-), Of, and r/. Hence, the network utility maximi- 
sation problem PI is not in the standard convex optimisation 
framework. Instead, we pose the following problem, 
P2: 



The sub-problem is given by 

rnax XI ^^ (^ ~ ^fi^f^ %' ^/)) 



subject to 



E </>/(A) = 1, V/G^ 

A=0 

0/(A) > 0, V/G J',VA<L>/-1. 



max max max V^ In (1 — ef{(f)f, Of,rf)) 

S> ft r ' ^ 



^ 9 r 



(4) 



/e^ 



subject to y^ 



kf 



f:ceCf 



'^fwf,c 



<Tc. 



\/ceC 



0/(A)>O, 
rf>Xf 



This is a separable convex optimisation problem, and hence 
can be solved by Lagrangian method. Let fif be a Lagrangian 

Df-l 

multiplier for the constraint ^ 0/(A) = 1, and define /j, = 

A=0 

[/i/]/^jr. The Lagrangian function is given by 
L($,A^) = Eln(l-e^)-E/iWl- E </>^(A) 

f ^ TT f^TT \ A— n / 



v/ G J^, < A < i:>/ 

V/G^ 
V/G^ 
V/G^ 

In general, the solution to P2 need not be the solution to 
PI. However, in our problem, we show that P2 achieves the 
solution of PI. 

Lemma 2. . For a function f : y x Z ^ M that is concave 
in y and in z, but not jointly in {y^ z), the solution to the joint 
optimisation problem for convex sets y and Z 



Applying KKT condition, 
1 



/e-F 



A=0 



is the same as 



max f{y,z) 

yey,zez 



(5) 



(6) 



maxmax/f?/, z), 

zez yey 

if f{y*{z)^ z) is a concave function of z, where for each z ^ Z, 
y*{z) :=argmax/(?/,2;). 

We note that for each r/ and Of, the probability of error 
Cf is convex in 0/, and hence, ln(l — e/) is concave in 4>f. 
Thus, we first solve for the optimum code bit allocation (j)% 
in Section IIV-AI Then, using the optimum code bit allocation, 
we solve for the optimum Chernoff bound parameter 6* which 
we describe in subsection IIV-BI After having solved for the 
optimum 6*, we show in Section HV-Cl that U{^*,6*{r), r) is 
a concave function of r. Hence, from Lemma [2l the solution 
to problem (P2) (the maximisation problem that separately 
obtains the optimum 0* and optimum r*) is globally optimum. 
We study the rate optimisation problem that obtains r* in 
subsection IIV-DI 



IV. Utility Optimum Rate Allocation 

A. Optimal Code Bit Allocation <l> 

We consider the maximisation problem defined in Eqn.|4]for 
a given coding rate vector r and Chernoff-bound parameter 
vector 6, and obtain the optimum 0/ for each flow f e T. 



dL 



d(l)f{i) 



U/(o* 



we get 



-e/ 



pfOfe^f^^f^'^ 



or, e 



Of^fii) 






/^/ 



(7) 



f3f i9/e/-/i/(l-e/) 

for i = 0, 1, 2, • • • , n/ — 1. Since, the RHS of Eqn. [7] is the 
same for all i, we get (/>j(^) = (t)%{j), and hence 

(/)}(A) = -^, VA = 0,1,---,I^/-L 

Thus, $* allocates equal portions of an encoded packet across 
transmission schedules with a delay of 0, 1, • • • , Df — 1, unlike 
the conventional transmission scheme which transmits all the 
coded bits of a packet in one shot. Hence, 6/(0^, Of^Vf) is 



Cf = exp 



0f{l-rf)-Df\n{l-l3f^l3fe- 



(8) 



B. Optimal 6>* 

We now consider the optimum Chernoff-bound parameter 
problem with the optimum coded bits allocation ^*, and for 
any given coding rate vector re [Aj, A/]^. 

max^ln(l-e^(0*,^^,r^)) (9) 

feT 
subject to <9/ > 0, V/ G J" 

We note that the objective function is separable in OfS, and 
that Cf is convex in Of. Hence, the problem defined in Eqn. ([9]), 
is a concave maximisation problem. The partial derivative of 
Cf with respect to Of is given by 



dcf 
80] 



-Cf 



(1 - rf) 



Pfe' 



f/Dj 



l-pf + jifc^i'^f 



Observe that 



Thus, if, for Of =0,1 



^^ — 0-jDj is an increasing function of Of. 



^f 



11. 



i--^7+^ < or r/ > 1 - /3/, 
the derivative is positive for all Of > 0, or e/ is an increasing 
function of Of. Hence, for r/ > 1 — (3f, the optimum 0*r is 
arbitrarily close to which yields e/ arbitrarily close to 1. 
Thus, for error recovery, for any end-to-end error probability 
f3f, the coding rate should be smaller than 1 — /3/, in which 
case, we obtain the optimum 6>j by equating the partial 
derivative of e/ with respect to Of to zero. 
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Thus, the probabiUty of a packet decoding error for a given 
r/ with the optimum allocation of coded bits #*, and the 
optimum Chemoff-bound parameter 9p is 



exp 



-Df 



(l-r/)ln 



rf 



Pf 



rf In 



rf 



1-/5/ 



exp{-DfKL{B{l-rf)\\B{l3f)) 



where KL(/i, /2) is the Kullback-Leibler divergence between 
the probability mass functions (pmfs) /i and /2. 

C. A convex optimisation framework to obtain optimal r% 

If ln(l — ef{4>*^^0*^^rf)) is concave in r/, then one can 
obtain the optimum r% using convex optimisation framework. 
To show the concavity of ln(l — e/(0j, ^p r/)) it is sufficient 
to show that ef{(f)pOp rf) is convex in r/. Note that 

6/ is convex if 

'1-r/ 



In 



Pf 



e/ 



In 



^/(l-^/) 



^/ 



1-/5/ 



> 



I^/ 



^/(l-^/)' 



or. 



/Df 



or, 



In 



In 



rf Pf 

l-Tfl-pf 



> 



/Df 



^rf(l-rf) 



<0 



^rf{l-rf) V rf (3f 

The function > ^ is convex in rf. Also, In ( ^^^ ) 

is decreasing with r/ , and hence, — In ( ^^ ~^ j 

/ l-A^ 1-/3 A 



In I 



Thus, we have a sufficient condition 



a/^/(i-^/) 



In 



1 - A/ 1 - /3/ 



A. 



<0 



(10) 



The above condition requires the delay deadline Df to be 
smaller than some Df{rf). We consider DfS to satisfy this 



condition, and hence, the rate optimisation problem is a 
concave maximisation problem. For the sake of completeness, 
we include this as a constraint in the problem formulation. 
However, this condition is not an active constraint. 



D. Optimal Coding Rate r 

From the previous subsection, we observe under the delay 
constraint Eqn. (fTOl) that ef{4>%^0*Arf)^rf) is convex in r/, 
and hence, we obtain the optimum coding rate r j using convex 
optimisation method. Also, from Lemma [2l it is clear that 
is the unique globally optimum rate. Thus, we solve the 



'/ 

following network utility maximisation problem 






\n{l-ef{ct>),0}{rf),rf)) 



(11) 



subject to 2_\ 



kf 



TfWf^, 



<Tr 



f:ceCf 

rf < A/ 
rf>Xf 

VD~f 



V^/_(l - rf) 
where a = In ( ^ ^ ~-^ 



a < 



yceC 

V/ G ^ (12) 



- o I . It is clear that the objective func- 
tion is separable and concave, and hence, can be solved using 
Lagrangian relaxation method. Also, we note here that the 
constraint represented by Eqn. ([T2l) is not an active constraint, 
and hence, there is no Lagrangian cost to this constraint. We 
note here that the coding rate should be such that kf /vf is 
an integer, and hence, obtaining r^ is a discrete optimisation 
problem. This is, in general, an NP hard problem. Hence, we 
relax this constraint, and allow Vf to take any real value in 
[A J, A/]. The Lagrangian function for the rate optimisation 
problem is thus 



L{r,p,u,v) 
^ln(l-e,)-^pj^ 



^/ 



VfWf^c 



feT cec \feTc 

+ Z! ^/ {rf - A/) - XI ^f i^f - ^/) 
feT feT 

Applying KKT condition, ^ |r*= 0, we have 
-1 dcf , v^ —kf 



ceCf 



l-cfdrf ■ ^ ^- rywf 



#Ei^ +»/-«/ 



/ \ceCf 



^/, 



If the optimum r% is either Aj or A/, then it is unique. If r^ G 
(A J, A/), then Uf = Vf = 0, which is the most interesting 
case, and we consider only this case for the rest of the paper. 



Let \f := E 



Pc 



. The above equation becomes 



6/ ,0* ^ ^f^f 



1-ef 



'f - ~~^ 



(13) 



e/ = T-r^^S^ (14) 



Xfkf + ey} 



* ^*2 



DfD{B{l-r})\\B{Pf))=\n 






(15) 



In the above equation, the LHS is a strictly convex decreasing 
function of r^. Since, the utiHty maximisation problem is a 
concave maximisation problem, the optimum rate r^ G (0,1 — 
Pf) exists and is unique. 



E. Sub-gradient Approach to Compute optimum p* 

In this section, we discuss the procedure to obtain the 
Shadow costs or the Lagrange variables p* . The dual problem 
for the primal problem defined in Eqn. ([TT1) is given by 

min D(p), 

p>0 

where the dual function D{p) is given by 
Dip) 

= max J2 ln(l " e/(r/)) + Y.P'^ [Te-Yl ^ ] 



(16) 



= ^ ln(l - ef{r}{p))) + ^p J Te - ^ 



kf 



/eJ^- 



cec 



feTa 



'*f{p)wf,c 



(17) 



In the above equation, ef{rf) denotes ef{<p*f,9*f{rf),rf). 
Since the dual function (of a primal problem) is convex, D is 
convex in p. Hence, we use a sub-gradient method to obtain 
the optimum p*. From Eqn. (fT6b . it is clear that for any r, 

D{p) > ^ ln(l - Cfivf)) + E^n ^^ - E r^ 



feJ^ 



cec 



feJ^c 



VfWf^c 



and in particular, D{p) is greater than that for r = rUp), i.e.. 
Dip) 

cec \ feTr. f^^^ -^'"^ 



feJ' cec \ feJ^c 

^(p) + E(^^-^^)(^^- E 



cec 



kf 
_^r}{p)wf. 



(18) 



Thus, a sub-gradient of £>(•) at any p is given by the vector 



Tr 



2-^ r* 



iP)Wf,c 



(19) 



cec 



We obtain an iterative algorithm based on sub-gradient method 
that yields p*, with p{i) being the Lagrangians at the ith 
iteration. 

+ 



Pc{i + 1) 



Pc{i) - 7 • M^c 






where 7 > is a sufficiently small stepsize, and [/(x)]+ := 
max{/(x), 0} ensures that the Lagrange multiplier never goes 
negative. Note that the Lagrangian updates can be locally done, 
as each cell c is required to know only the rates r*^{p{i)) 
of flows f ^ Tc- Thus, at the beginning of each iteration z, 
the flows choose their coding rates to rj(p(i)), and each cell 
computes its cost based on the rates of flows through it. The 
updated costs along the route of each flow are then fed back 
to the source node to compute the rate for the next iteration. 
The Lagrange multiplier p^ can be viewed as the cost of 
transmitting traffic through cell c. The amount of service 
time that is available is given by (5 = Tc — ^ 



feT, 



r*(p(i))w/,, 



When b is positive and large, then the Lagrangian cost pc 
decreases rapidly (because D is convex), and when (5 is 
negative, then the Lagrangian cost Pc increases rapidly to make 
J > 0. We note that the increase or decrease of pc between 
successive iterations is proportional to (5, the amount of service 
time available. Thus, the sub-gradient procedure provides a 
dynamic control scheme to balance the network loads. 

We explore the properties of the optimum rate parameter 
r J in Section IIV-F[ In Section |Vl we provide some examples 
that illustrate the optimum utility-fair resource allocation. 

F. Properties of t% 

Lemma 3. t%(Dj) is an increasing function of Df. 

Lemma [3] is quite intuitive. For any given channel error /3/, 
as the deadline become less stringent, it is optimal to go for a 
high rate code. In other words, it is optimal for a flow to use as 
much scheduling time as possible (for a large Df, and hence, 
use a high rate code); however, the resources are shared among 
multiple flows, and hence, we ask the following question: 
''what is the optimal share of the scheduling time'' that each 
flow should have. Interestingly, in our problem formulation, 
the code rate r/ also solves this optimal scheduling times for 
each flows. 

V. Examples 

A. Example 1: Two cells with equal traffic load 

We begin by considering the example shown in Figure H] 
consisting of two cells a and h having three nodes 1, 2, and 
3. Each cell has the same packet erasure probability P and the 
schedule length T. There are three flows /i, /2, and /s, with 
two of the flows /i and /a having one-hop routes C/^ = {^} 
and C/3 = {a}, and one flow /2 having a two-hop route 





Fig. 4. Cells with equal traffic load 



Fig. 5. Cells with unequal traffic load 



C/2 = {a, 6}. Each flow has the same information packet size 
k, decoding deadhne D and PHY transmit rate, i.e. Wf^c = ^• 
This is analogous to the so-caUed parking-lot topology often 
used to explore fairness issues. 

The end-to-end erasure probability experienced by the two- 
hop flow /2 is greater than that experienced by the one hop 
flows /i and /a, since each hop has the same fixed erasure 
probability. Hence, we need to assign a lesser coding rate rf^ 
to flow /2 than to flows /i and /s in order to obtain the 
same error probability (after decoding) across flows. However, 
when operating at the boundary of the network capacity region 
(thereby maximising throughput), decreasing the coding rate 
r/2 of the two-hop flow /2 requires that the coding rate of 
both one-hop flows /i and fs be increased in order to remain 
within the available network capacity. In this sense, allocating 
coding rate to the two-hop flow /2 imposes a greater marginal 
cost on the network (in terms of the sum-utility) than the one- 
hop flows, and we expect that a fair allocation will therefore 
assign higher coding rate to the two-hop flow /2 . The solution 
optimising this trade-off in a proportional fair manner can be 
understood using the analysis in the previous section. 

In this example, both the cells are equally loaded and, by 
symmetry, the Lagrange multipliers pa = Pb- Hence, A/^ = 
-^ = A/3. For the Chernoff-bound parameter 6 = [0^0], we 
find from Eqn. ([13]) , 
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For sufficiently small erasure probabilities, we have 
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flows while cell a carries 


only one flow. The encoding rate 


constraints are given by 
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Since, both r/^ and r/2 are at most 1, it is clear that at 
the optimal point, the rate constraint of cell a is not tight 
while the constraint of cell b is tight. Thus, the shadow prices 
(Lagrange multipliers) Pa = and ^5 > 0. That is, at the first 
hop the cell is not operating at capacity, and so the "price" 
for using this cell is zero. In this example, A/^ = A/2, and 
hence, from Eqn. ([T3l) . we deduce that for sufficiently low 
cell erasure probability /3, e/^ ^ Cf^. Alternatively, as the 
delay deadline D ^ 00, from Eqn. iTSh we have e/^ = Cf^. 
These proportional fair allocations make sense intuitively since 
although flow /2 crosses two hops, it is only constrained at the 
second hop and so it is natural to share the available capacity 
of this second hop approximately equally between the flows. 
When the erasure probability is sufficiently small, this yields 
approximately the same error probabilities for both flows. For 
larger erasure probabilities, it leads to the two-hop flow having 
higher error probability, in proportion to the per-hop erasure 
probability (3. 
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Thus the proportional fair allocation is e/^ = e/3 ^ 1/2 • Cf^. 
That is, the coding rates are allocated such that the one-hop 
flows have approximately half the error probability of the two- 
hop flow. 



B. Example 2: Two cells with unequal traffic load 

We consider the same network as in the previous example, 
but now with only the flows /i and /2 (i.e., the flow fs is not 
present) in the network. In this example, cell b carries two 



